---
abstractspacing: double
fontsize: 12pt
margin: 2cm
urlcolor: darkblue
linkcolor: Mahogany
citecolor: Mahogany
spacing: single
bibliography: references.bib
biblio-style: apalike
output:
  pdf_document:
    citation_package: natbib
    fig_caption: no
    number_sections: no
    keep_tex: no
    toc: no
    toc_depth: 3
    template: article-template.latex
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, warning=FALSE, message=FALSE,cache=TRUE)
knitr::opts_chunk$set(fig.width=7, fig.height=6, out.width = '70%', fig.align = "center") 
rm(list=ls())
library(remotes)
library(kableExtra)
library(haven)
library(tidyverse)
library(ivmodel)
library(doParallel)
library(foreach)
library(estimatr)
require(AER)
library(lfe)
library(glue)
#path <- "/Users/ziwenzu/Dropbox/research/IV/IV Sensitivity/LLXZ_rep"
path <- "~/Dropbox/ProjectZ/IV Sensitivity/LLXZ_rep"
setwd(path)
knitr::opts_knit$set(root.dir = path)
#install_github("apoorvalal/ivDiag")
library(ivDiag)
# number of cores
cores <- 15
```


## @tajima2013

| Replication Summary | |
|----------------|----------------------|
| Unit of analysis | village and urban neighborhood |
| Treatment | distance to police posts (as a proxy for exposure to military intervention) |
| Instrument | distance to health station |
| Outcome | incidence of communal violence |
| Model | Table1(4)|

```{r ajps_Tajima_2013}
df<-readRDS("./rawdata/ajps_Tajima_2013.rds")
D <-"z2_distpospol"
Y <- "horiz2"
Z <- "z2_dispuskes"
controls <- c("flat", "z2_altitude","urban", "natres", "z2_logvillpop", "z2_logdensvil",
              "z2_povrateksvil", "z2_fgtksvild", "z2_covyredvil", "z2_npwperhh", 
              "z2_ethfractvil","z2_ethfractsd", "z2_ethfractd", "z2_relfractvil", 
              "z2_relfractsd", "z2_relfractd", "z2_ethclustsd", "z2_ethclustvd", 
              "z2_relclustsd", "z2_relclustvd", "z2_wgcovegvil", "z2_wgcovegsd", 
              "z2_wgcovegd", "z2_wgcovrgvil", "z2_wgcovrgsd", "z2_wgcovrgd",
              "natdis","javanese_off_java", "islam", "split_kab03", "split_vil03")
cl <- 'kabid03'
FE <- 'prop'
weights<-"probit_touse_wts03"
(g<-ivDiag(data=df, Y=Y, D=D, Z=Z, controls=controls, FE =FE, 
  cl =cl,weights=weights, cores = cores))
```

```{r, cache = FALSE}
plot_coef(g)
```

```{r ajps_Tajima_2013_sav, echo = FALSE}
save(g, file="./estimate/Tajima2013.RData")
```